Comment on "On the ionization equilibrium of hot hydrogen plasma 
and thermodynamic consistency of formulating finite partition functions" 
[Phys. Plasmas 17, 062701 (2010)] 
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(Dated: March 8, 2011) 

Zaghloul [Phys. Plasmas 17, 062701 (2010)] reconsiders the occupation probability formalism in plasma thermody- 
namics and claims inconsistencies in previous models. I show that the origin of this incorrect claim is an omission of 
the configurational factor from the partition function. This arXiv version is supplemented with two appendices, where 
I add remarks and comments on two more recent publications of the same author on the same subject: on his response 
to this Comment [Phys. Plasmas 17, 124705 (2010)] and on his criticism towards the Hummer and Mihalas's (1988) 
formalism [Phys. Plasmas 17, 122903 (2010)]. 

PACS numbers: 52.25.Kn, 05.70.Ce 

In a recent paper, Zaghloul 1 revised the occupation prob- 
ability formalism routinely applied for quenching divergen- 
cies in frames of the chemical picture of plasmas. 2,3 Following 
Ref.[H he considers a plasma composed of protons, electrons, 
and H atoms and writes separate expressions for the contribu- 
tions of these subsystems into the free energy: F e , F p , and Fa, 
respectively. The atomic contribution is written in the form 



Fr = N H k B T 



In 



(1) 



where k B is Boltzmann constant, T temperature, Nr the 
total number of atoms in all quantum states, Ah = 
(2ttH 2 /mk B T) 1 / 2 the thermal wavelength of an atom, m is 
the atomic mass, and Qi n t,H is the internal partition function. 
The author fails to notice that Eq. ([l} is valid only for a Boltz- 
mann gas of noninteracting particles (e.g., Ref. 4, §§41, 42). 
In general, instead of Eq. ([TJ one should start from the ex- 
pression F = —k B T Tr e~ H , where H is the total Hamilto- 
nian of the system (e.g., Ref.0, § 31). Assuming that (i) the 
motion of particles is quasi-classical, (ii) the kinetic and po- 
tential energies in H are uncoupled, (iii) interactions between 
plasma particles appear in H as an additive potential function, 
one has 5 ' 6 F = -k B TlnZ = -k B Tln(Z tlans Z int Z con{ ) = 
-Ftrans + -Pint + F con f, where the first two terms correspond 
to the translational and internal degrees of freedom and the 
third one takes into account interactions between all plasma 
particles (in general, not only those between neutral atoms). 
In the case of H atoms, lnZ trans , H = -F tlans . B /k B T = 
NjihiieV/NjiX^). Having defined Q conf = Z { 
Qint = Zfa^" , one can write 



1/JVh 



and 



F H = N H k B T 



In 



N B X 3 H 



VQ int Q 
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(2) 



In general, Eq. © cannot be reduced to Eq. ([T). Moreover, 
since level populations depend on interactions in the plasma, 



<2i„t in Eq. (f2| may differ from Q- lrA h for the ideal Boltz- 
mann gas in Eq. (Q~|) (it is well known 2 - 5 that Qi n t,H needs a 
cutoff to avoid divergency due to the infinite number of shal- 
low Rydberg states). Conversely, <5 con f depends on internal 
level populations, because interaction forces between atoms 
depend on their excitation states. Thus, Fa& and F con { are not 
independent, and the definition of F- nlt is not obvious. 

The free energy minimization method assumes that F is 
expressed explicitly through numbers of particles of differ- 
ent kinds and minimized with respect to these numbers at 
constant volume V. In our case, F = F({N K }, N e , N p ), 
where A^ are numbers of atoms on quantum levels k. Let 
us calculate Fid = -Ft rans + Fj nt using relation 4 F = 
E — TS, where E is the mean energy and S is the en- 
tropy. Assuming that the plasma is uniform in space, and 
motion of atoms is classical with distribution density T K (p) 
over momenta p, the contribution of A^ K atoms to E is 
A'k J d 3 p J~k{p) e«(p)) where e K (p) is the total (kinetic mi- 
nus binding) atomic energy, while the entropy contribution is 
-k B N K Jd 3 P F K (p) ln[F K (p) (2nh) 3 Njg K eV)} where g K 
is quantum degeneracy of level k. Let us consider the case 
where e K (p) = p 2 /2m — \k and binding energies \k do 
not depend on p (a more general case has been studied in 
Ref.0). ThenJ" K (p) = (A H /27rft) 3 e" p2 / 2m ' £BT . After in- 
tegration and adding the translational contribution of N p clas- 
sical protons and the contribution of electron gas i*id, e , one 
obtains 



F id = k^Y^N^He-^I^^N^l/g^V) 

K 

+k B TN p [ln(N p X 3 p /V) - l] + F idj0 , (3) 



where X p is the proton thermal wavelength. For brevity we 
shall approximate A p = Ah- The minimum of F = Fid + 
i^onf under the stoichiometric constraints with respect to dis- 
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sociation/recombination reactions H ^ e + p requires 

dF dF dF 



dN K dN p dN e 

This gives, with account of Eq. (01, 



In 



N K /g K 



k B T 



df , df df 



dN p dN e dN K ' 



(4) 



(5) 



where fi e = dF id , e /dN e and / = F con{ /k B T. 

An occupation probability w K is conventionally defined 2 as 
the probability of finding the atom in state k relative to finding 
it in a similar ensemble of noninteracting ions. In our case this 
means that N K oc w K g K e XK ^ kBT . Therefore, according to Eq. 
(01, In w K = —df/dN K + Ch, where Ch does not depend on 
N K . Thus one can write 



N K _ w K g K e 
N B ~ 



where 



Qi 



int.H,' 



int, YL.w 



(6) 



(7) 



Note that number fractions N K /N B do not depend on Ch- 
Hummer & Mihalas 2 set Ch = 0. However, an additional re- 
quirement that the equation of ionization equilibrium for non- 
degenerate plasma has the form of Saha equation multiplied 
by w K [N K oc N p N e w K e^l k * T ; seeEq. (17) of Ref.H leads 
to 



In w K 



df , df df 



8N n 8N B dN K 



Ch,, 



(8) 



where Ch e p is independent of N K , N e , and N p . Given the 
constraints N B = ^ K an d + N p =constant, it is easy 
to see that N K do not depend on the choice of C B , e ,p- We set 3,7 
C H ,e, P = (then obviously C H = df/dN p + df'/dN e ). 
Substitution of (01 into (0 gives 

k B TN B [ln{N B \yV) - 1] 
+k B TN p []n(N p X 3 p /V) - 1] + F id , c + F lnt , (9) 



F 



id 



where 



Fmt = -k B TN B In Qint.H.w + k B T ^ N * m 



(10) 



Note that Qi n t t n >w appears in (O merely as a normaliza- 
tion constant, and the occupation probabilities w K are auxil- 
iary quantities, defined from the condition of the minimum of 
the total free energy according to Eq. (|8)- 

Zaghloul 1 follows another route. He replaces Qint.H by 
Qint,H,ii> m Eq- leaving the meaning of quantities w K un- 
defined, and assumes that this replacement is a way of ac- 
counting for the nonideality effects, alternative to the intro- 
duction of F con f (as he explicitly writes and exposes in his 
Eq. 26). This implies that the product Qi n tQ C onf in Eq. (O 
can be represented as a single sum In general, it cannot. 



Furthermore, this assumption leads to an additional restriction 
on w K (Eq. 32 of Ref. 1), which may not necessarily be ful- 
filled in a real plasma. 

We should remark that the expression for the free en- 
ergy can be written through w K without F con f in the "low- 
excitation approximation" of Hummer & Mihalas, 2 who write 
it in the form / — ^ K N K df /dN K = 0. Taking into account 
that they consider the case where Ch = 0, this approximation 
can also be written as 



F 



conf 



k B Tj2 N * 



In ) 



= 0. 



(11) 



The latter form is more general. When condition (ITT1 is 
satisfied, the second term in Eq. ( TTOb annihilates with the 
configurational part F con f of the total Helmholtz free energy 

F ^trans ~t~ ^int ~t~ -^conf ■ 

The low-excitation approximation has serious shortcom- 
ings (see discussion in § Hid of Ref. 0). One can explic- 
itly show that it is violated in some thermodynamic models 
commonly used in literature (for instance, the hard-sphere 
model 2 ). For these reasons, approximation ([Til is used rather 
rarely. In particular, it was not employed in Refs. H] and 
Without this approximation, however, F = F^ + F con { does 
not reduce to an expression containing only w K without F con f , 
as required in Ref. |lj 

In short, the conclusions in Ref. 1 originate from a trivial er- 
ror: the author arbitrarily removes from the partition func tion 
the configurational factor that is responsible for interactions 
between plasma particles, however assumes the significance 
of such interactions by allowing occupation probabilities to 
differ from unity. The controversies in Ref. [l] result from this 
basic omission and not from the alleged inconsistencies of the 
previous models. 

This work was partially supported by Rosnauka Grant NSh- 
3769.2010.2 and RFBR Grant 08-02-00837. 



Appendix A: Comment on Zaghloul's response [Phys. 
Plasmas 1 7, 1 24705 (201 0)] 

After the publication of this Comment 8 , Zaghloul 
responded 9 with twelve items of arguments, one of which is 
subdivided into four subitems. Here I briefly comment on 
them, item by item. 

(i) "Equation (2) can always be reduced to Eq. (1) by defin- 
ing QintQconf = Qint,H,u) " — but in this case Q int ,H,w will 
not have the required form of a partition function with weights 
w K , Eq. (0. 

(ii) " Following Eq. (2), the author writes: 'Moreover, since 
level populations depend on interactions in the plasma, Q- m t 
in Eq. (0 may differ from Qi n t, B for the ideal Boltzmann gas 
in Eq. (0 (it is well known^ that Qi n t,H needs a cutoff to 
avoid divergency due to the infinite number of shallow Ryd- 
berg states). Conversely, Q CO nf depends on internal level pop- 
ulations, because interaction forces between atoms depend on 
their excitation states. Thus, F ln t and F con f are not indepen- 
dent, and the definition of F mt is not obvious.' Clearly, this 
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is the main message of our paper." — Here Zaghloul agrees 
that i*i nt and F con f are not independent. However, next he 
writes: "according to factorizability of the PF adopted in [the 
Comment], Q m t in Eq. (2) should be Qi n t,H as in Eq. (1)." — 
This is not logical. Just the opposite, since Fj nt and F con { are 
not independent, the presence of nonideality (i.e., F con f ) mod- 
ifies level populations and thus stipulates the modification of 
<9int,H which transforms into Qi n t,H,w 

(iii) Having cited Eq. 0, Zaghloul then writes: "the 
above definition of w K can be written as w K — 
{N^ W /N R )/(N K /N H ) = N K , W /N K ." — This equation is 
mathematically wrong, as its left-hand side does not agree 
with Eq. 0. Therefore, the subsequent discussion in this item 
of is incorrect as well. 

(iv) "The substitution of Eq. (6) into Eq. (3) [...] repre- 
sents an inconsistency [...], where Eq. (3), which embodies 
uncoupling between the translational and internal energies, 
is derived using a Boltzmann distribution of the excited states 
while Eq. (6) represents a different distribution (real distribu- 
tion [...])" — in fact, however, Eq.(3) is derived using the 
real distribution (which has to be determined); it is not the 
ideal-gas Boltzmann distribution. 

(v) Discussing the last term on the right-hand side of Eq. 
( II Oi l. Zaghloul writes: "The negative of this term, namely, 
~^bT2~2 k N k In w K , was interpreted [. . . ] as a contribution 
to the ideal-gas part of the entropy due to correction w K to 
the probability that the nth state is occupied. However, this 
is physically incorrect since w K is less than unity by defini- 
tion. As a result, the contribution stated above is positive 
which means [. . . ] that using the occupation probabilities 
w K , which diminishes levels' degeneracies, will increase the 
entropy, which is in direct contradiction with the fundamental 
statistical interpretation of the entropy!." — First, in frames 
of the formalism under consideration 2 ' 3 (where, for simplic- 
ity, quantum-mechanical bound level shifts and broadening in 
the plasma environment are neglected), using the occupation 
probabilities w K does not "diminish the levels' degeneracies": 
quantum-mechanical degeneracies remain the same, g K . In 
this approach, the occupation probabilities are additional fac- 
tors, which have a different physical meaning. Second, the 
factor w K , as defined in Eq. 0, is the ratio of nonideal to 
ideal occupation numbers, which is not necessarily less than 
unity (see discussion in Ref. [It). Third, the increase of the 
entropy with deviation of the true distribution of level popu- 
lations from the Boltzmann distribution is not in "contradic- 
tion," but in full accord with fundamental principles of statisti- 
cal thermodynamics. Indeed, since the Boltzmann distribution 
provides the minimum entropy (without interactions), then the 
use of w K must increase the part of the entropy related to Fi nt . 

(vi) This item attempts to disprove Eq. and consists of 
four subitems ( "several mistakes in the derivation of this ex- 
pression need to be unveiled here", writes Zaghloul). They 
are commented in order below. 

(vi-1) Equation is "based on the above-quoted defini- 
tion of w K . We have shown above that this definition is inac- 
curate." — Naturally, the equation is based on the definition, 
there is no contradiction. As concerns the alleged inaccuracy 
of this definition "shown above," this statement arises from 



Zaghloul's confusions listed (as well) above. 

(vi-2) Equation "is in every respect equivalent to Eq. 
0. Since no efforts were devoted in the derivation ofEq. 
to assure convergence of the IPF, then there is no guarantee 
for the convergence of the IPF" — this statement is true, but 
there is no contradiction, nor inconsistency. Equations and 
are quite general, while the convergence depends on the 
existence of i^onf and on its form. 

(vi-3) "no efforts were made to assure that N K /Nn satisfies 
the normalization condition [Eq. 07" — this point is simply 
wrong, because, given Eq. 0, the normalization condition 
2^2 K N K = iVn is equivalent to definition 0. 

(vi-4) "Most importantly, we have shown mathematically 
[. . . ] that for occupational probabilities similar to those given 
by Eq. [. . . ] the condition for the separability of the con- 
figurational component [.. . ] implies a linear dependence of 
F CO ni on the populations of individual excited states " — in 
fact, this "mathematical" demonstration was done in frames 
of the assumption that F con f should be dropped out as soon as 
Fi nt is modified. The failure of this assumption is just in focus 
of my Comment. 

(vii) This item consists of three paragraphs. The first one 
contains two long quotations, one from the Comment and an- 
other from the work by Hummer & Mihalas 2 , both concerning 
the low-excitation approximation. It is followed by two para- 
graphs of the criticism of the Hummer & Mihalas's work with 
a reference to another paper by Zaghloul 10 , especially devoted 
to this criticism. The incorrectness of the latter paper 10 is ex- 
plained in AppendixlBlbelow. 

(viii) This item consists of three paragraphs filled with ex- 
cerpts from the last part of the Comment, which deals specifi- 
cally with the low-excitation approximation. Zaghloul claims 
it to be "completely irrelevant as we do not use such an ap- 
proximation, and in fact we strongly criticized it", but fails to 
note that his assumption of the cancellation of F con f with the 
last term in Eq. (fTOT t is equivalent to Eq. (fTTT i. that is equiv- 
alent to using the same low-excitation approximation that he 
pretends to criticize. 

The last paragraph of item (viii) ends up with the follow- 
ing passage: "Axiomatically, the general case embodies the 
special case as one of its possibilities. Now, is it possible that 
the special case shows a restriction that the general case does 
not show?" — Here, the special case is meant to be the low- 
excitation approximation. Naturally, the special case always is 
a restricted applicability relative to the general case, therefore 
this rhetoric question is pointless. 

(ix) This item consists in a reformulation of Eq. (though 
with some misprints), followed by the sentence: "It is left to 
the reader to verify that the condition of the minimum" [of the 
free energy, Eq. 0] "cannot lead to the expression for w K 
given by Eq. 0." — In fact, the reader can easily see that Eq. 
is the direct consequence of the condition of the minimum 
of the free energy, as Eq. is derived from Eq. and Eq. 
0, through Eq. 0, using definitions and 0. 

(x) The author refers to his Ref. [l0| (to be dealt with in Ap- 
pendixIBl. quotes the expression w = exp(— F con { / N k^T) , 
and says: "the claim" [in the Comment] "that an error origi- 
nates in our analysis from removing the configurational com- 



4 



ponent from the PF and assuming the significance of such in- 
teractions by allowing occupation probabilities to differ from 
unity, is groundless." — In fact, this sentence is groundless it- 
self, because the Comment shows that Zaghloul's omission of 
-Pconf from the final expression of F (after Fj nt has been mod- 
ified by introducing w K ) is incorrect, in general, unless one 
additionally uses the low-excitation approximation, Eq. (fTTT >. 

(xi) The author states that summing over k in Eq. <[8J 
amounts to "the term commonly referred to in the literature as 
a lowering of ionization energy" — which only confirms the 
consistency the results exposed above with other well known 
results in the literature. Zaghloul writes further: "Clearly, the 
inclusion of this term cannot exclusively lead to a truncation 
of the IPF", which is a mistake: it can lead to such truncation. 

(xii) "For completeness, we would like to explain here that 
Eq. (0, or more generally the relation Fn = —ksThiZ is 
valid for Boltzmann distribution for the total energy [.. .] of 
the system " — The latter relation, where Z is the total parti- 
tion function of the system, including all ideal and nonideal 
contributions, is a textbook fact. However, the general rela- 
tion Fh = — ksT In Z is not equivalent to its particular case, 
Eq. (fl}. In the Comment I called to the reader's attention, that 
one should not substitute into this relation Z = Z tTans Zi nt , 
omitting Z coa {, as Zaghloul does. 

Appendix B: Comment on another similar paper by 
Zaghloul [Phys. Plasmas 17, 122903 (2010)] 

A more recent paper 10 of the same author is devoted to the 
criticism of the formalism of Hummer & Mihalas^, which lies 
in the basis of the Opacity Project (e.g., Ref. Ull an d refer- 
ences therein). There are several mistakes in Ref. 10. One of 
them is that the author mixes up the occupation probabilities 
w K with quantum statistical weights of the bound states g K , 
writing Wi — giUJi = exp[— {df /dNi)/k^,T\. The statistical 
weight gi arises from counting multiply degenerate quantum 
states, gi ^ 1, whereas exp[— (df /dNi)/k-QT] is most of- 
ten less than 1. This mistake has no effect, for example, if one 
considers the case of = 1. Another mistake is the statement 
that "a monatomic perfect gas has only translational kinetic 
energy ( no means of storing energy except as kinetic energy)." 

However, there is a more serious fault that lies in the basis 



of the article and hence disproves its main conclusions. 

The main focus of the criticism in this paper is the expres- 
sion F = -Ft rans + F[ nt + F con {. The author writes: "the 
factorizability of the partition function ( or, equivalently, the 
separability of the free energy components) implies that vari- 
ous types of energies are independent of each other" . In fact, 
there is no such implication. He continues: "The separation 
of the configurational free energy therefore indicates that it 
has no influence on the internal free energy." In fact, it de- 
pends on the definition of the Fi nt and F con { terms, which is 
not obvious in a nonideal system. If F; nt is defined using the 
real occupation numbers (as done, e.g., in Refs. 3,7 and in 
the Comment above), then F con { affects the value of Fi nt at 
equilibrium through N K values. In general, the minimum of 
F = F int + F trans + F conf realizes at a different set of {N K } 
than the minimum of F nt + F tlans alone. In other words, the 
separability does not imply independence. 

Furthermore, the author writes: "if the configurational free 
energy ( or interaction energy) has no influence on the internal 
free energy, the expectation that including a separable config- 
urational component could lead to a truncation of the internal 
partition function [. . . ] is physically and logically question- 
able" (in the e-print version [arXiv:1010.1102vl], the word 
"questionable" is replaced by "incorrect"), "because they are 
independent of each other by assumption" . However, there is 
no such assumption in the cited references. On the contrary, 
since both F- mt and F con { are functions of the particle numbers 
(including occupation numbers), they are interrelated. 
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